Description of nuclear octupole and quadrupole deformation close to the axial 
symmetry and phase transitions in the octupole mode 
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The dynamics of nuclear collective motion is investigated in the case of reflection-asymmetric 
shapes. The model is based on a new parameterization of the octupole and quadrupole degrees of 
freedom, valid for nuclei close to the axial symmetry. Amplitudes of oscillation in other degrees 
of freedom different from the axial ones are assumed to be small, but not frozen to zero. The 
case of nuclei which already possess a permanent quadrupole deformation is discussed in some 
more detail and a simple solution is obtained at the critical point of the phase transition between 
harmonic octupole oscillation and a permanent asymmetric shape. The results are compared with 
experimental data of the Thorium isotopic chain. The isotope 226 Th is found to be close to the 
critical point. 
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I. INTRODUCTION 

Phase transitions in nuclear shapes have been recently 
observed at the boundary between regions characterized 
by different intrinsic shapes of quadrupole deformation. 
It has been shown by Iachello P, 0> tl that new dy- 
namic symmetries, called E(5), X(5) and Y(5) hold, re- 
spectively, at the critical point between spherical shape 
and 7-unstable deformation, between spherical and de- 
formed axially symmetric shape and between deformed 
axial and triaxial shape. Here, we are mainly inter- 
ested in the second case. First examples of X(5) sym- 
metry in transitional nuclei have been found in 152 Sm 
0| and 150 Nd Q. Other candidates for the X(5) sym- 
metry in different nuclear regions have been reported 

later HHHHHIUHIS. Fr0m the theoretical 
point of view, slightly different potentials have been ex- 
plored by Caprio | 14| and by Bonatsos et al. .131, while 
the evolution of the s = 2 band when the lower border of 
the square well potential is displaced from zero has been 
investigated by Pietralla and Gorbachenko 0] . 

A similar phase transition (i.e., from shape oscillation 
to permanent deformation, conserving the axial symme- 
try of the system) could take place also in the octupole 
degree of freedom. We have found a possible example 
of such a phase transition in the Thorium isotope chain, 
with the critical point close to the mass 226. In this case, 
the octupole mode is combined with a stable quadrupole 
deformation. Preliminary results have been reported in 
two recent Conferences 

sua. 

In order to provide a theoretical frame where to discuss 
the different aspects of the octupole motion, we introduce 
a new parameterization of the quadrupole and octupole 
degrees of freedom, valid in conditions close to the axial 
symmetry. In this limit, a model similar to the classical 
one by Bohr [l^ has been developed. This is the subject 
of the first part of this paper. In the second part, the 



model is used to discuss the evolution of the octupole 
mode along the isotopic chain of Thorium, and the results 
are compared with the experimental data. 



II. THEORETICAL FRAME FOR COMBINED 
QUADRUPOLE AND OCTUPOLE EXCITATIONS 

A. Previous investigations of the Octupole plus 
Quadrupole deformation 

Reflection-asymmetric nuclear shapes have been dis- 
cussed in a number of papers, either in terms of sur- 
face Quadrupole + Octupole deformation (Bohr geomet- 
rical approach) or with an extended Interacting Boson 
Model (algebraic approach) . The latter, proposed in 1985 
by Engel and Iachello [l!j, has been recently used by 
Alonso et al . I2f3 . Raduta et al. 0, l^, Zamfir and 
Kusnezov [2a , |24| . An alternative approach assuming a- 
cluster configurations has been discussed by Shneidman 
et al. 25]. In the frame of the geometrical approach, a 
number of theoretical investigations of the octupole vi- 
brations around a stable quadrupole deformation have 
been reported in the last 50 years |H EHH |H US El 
132] ]. Most of them, however, are limited to the case of 
axial symmetry. This ap pro ach has been criticized, e.g., 
by Donner and Greiner |33j |. who have stressed the fact 
that all terms of a given tensor order must be taken into 
account for a consistent treatment. To do this, Donner 
and Greiner renounce to the use of an "intrinsic frame" 
referred to the principal axes of the overall tensor of in- 
ertia and choose to define the octupole amplitudes in 
the "intrinsic frame" of the quadrupole mode alone. In 
this approach, definite predictions have been obtained at 
the limit where the octupole deformations are "small" in 
comparison with the quadrupole ones |34| . 
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B. The new parameterization 

Here we adopt a different approach, which can be 
useful also in the case of comparable octupole and 
quadrupole deformation, close to the axial-symmetry 
limit. Namely, we choose as "intrinsic" reference frame 
the principal axes of the overall tensor of inertia, as it 
results from the combined quadrupole and octupole de- 
formation. The definitions of quadrupole and octupole 



amplitudes ajti 



with A = 2,3 and = 
are recalled in the Appendix ^] All these amplitudes 
are defined in the (non inertial) intrinsic frame. To this 
purpose, in the case of the quadrupole mode alone, it is 



enough to assume — rea ^ an d 



„(2) 



: 0, with the 



standard parameterization in terms of 02 and 72: 

J 2 ) - 



,(2) 



(2) 




(la) 
(lb) 
(1c) 



For the octupole mode alone, a parameterization suitable 
to this pur pose has been proposed in 1999 by Wexler and 
Dussel |33 . We adopt here a very similar one: 



4 3) = 

4 3) = 

i (3) - 

,( 3 ) _ 



/3 3 COS 73 

-(5/2) (X + iY) sin 73 

y/l/2 (3 3 sin 73 

X COS73 + (VT5/2) sin73 

i Y cos 73 — (Vl5/2) sin 73 



(2a) 
(2b) 
(2c) 

(2d) 



With this choice, the tensor of inertia turns out to be 
diagonal (see Appendix lA"ll. 

In both cases, one has to consider, in addition to 
the intrinsic variables (P2, 72 for the quadrupole, or 
/?3 , 73 , X , Y for the octupole) , the three Euler angles 
defining the orientation of the intrinsic frame in the lab- 
oratory frame, in order to reach a number of parameters 
equal to the number of degrees of freedom (5 for the 
quadrupole, 7 for the octupole). 

Unfortunately, the situation is not so simple when 
quadrupole and octupole modes are considered together, 
as the intrinsic frames of the two modes do not necessar- 
ily coincide. We shall limit our discussion to situations 
close to the axial symmetry limit - in which, obviously, 
the two frames coincide - and define a parametrization 
which automatically sets to zero the three products of 
inertia 3 K , K < (k 7^ K,') up to the first order in the ampli- 
tudes of non-axial modes. 

To this purpose we put 



,(A) 



where are defined according to the eq.s 11121) and a 



4 A) 

are correction terms, which are assumed to be small com- 
pared to the axial amplitudes , but of the same order 



of magnitude as the other non-axial terms. It will be 
enough to consider these corrections only for those am- 
plitudes which, according to eq.s 11121) . arc cither zero or 

small of the second order: the imaginary part of a. 

and the real and imaginary parts of a x 



(2) (3) 



(2) 
2 > 



,(3) 



ft. 



The six 

"new" first-order amplitudes added to those of eq.s ijll^l) 
are, however, not independent from one another, if we 
choose as the reference systeme the one in which the three 
products of inertia turn out to be zero. 

The expressions of the inertia tensor as a function of 
the deformation parameters, obtained with the Bohr as- 
sumptions of not-too-big deformations and irrota- 
tional flow, are given in the Appendix ^ In order to 
simplify the notations, from now on we consider the in- 
ertia parameters B2 , -B3 included in our definitions of the 
amplitudes , which therefore correspond to \JB\ 
in the original Bohr notations. From the eq.s (|A20|I . and 
retaining only terms of the first order in the small ampli- 
tudes a}!/ with fi 0, we obtain the conditions 



i7l2 

J13 



-2V6 (/3 2 Im 



;(2) 



V5/3 3 Im 4 3) 



«J 2 3 



(foa? ) + V2/3 3 ai 3) 



: (4a) 
(4b) 



which are satisfied (at the leading order) if we put 



V2/3.3 



(v + <) ar = 



= (3) §2 



Im a 



~(2) 



V503 



Im a 2 3 ^ 



(v + <) 



(5) 



with the new parameters 77, £ and £ small of the first 
order. It is clear that only the ratios of the relevant am- 
plitudes are constrained by the eq.s (QJ. The definition 
of the new variables given in each line of eq.s Q con- 
tain therefore an arbitrary factor. Our choice (and in 
particular for the square-root factors at the denomina- 
tors) has some distinguished advantage, which will result 
clear from the classical expression of the kinetic energy, 
discussed in the next paragraph. 

In the intrinsic reference frame, and at the same or- 
der of approximation, the values of the three principal 
moments of inertia can be derived from eq.s l|A18IA19() : 



Ji 
Ji 
J?, 



2/91) 



2V3(/3| 7 2 



3(/? 2 2 



2/? 2 ) - 2V3(/3 2 2 72 - 
4(/?2 2 7 2 2 + /?3 2 7 3 2 ) + 18(A 2 
+2{ v 2 + C 2 ) + 8£ 2 . 



\/5/3 3 2 73) 

V5/3 3 2 7 3) 
Y 2 ) 



(6a) 
(6b) 

(6c) 



With the amplitudes given by eq.s 11121) the principal 
axes of the quadrupole would coincide with those of the 
octupole. It is not necessarily so with our more general 
(3) assumptions. When a[ X ^ 7^ 0, the axis 3 of the tensor of 



inertia for the quadrupole mode alone does not coincide 
with that of the octupole. If a\ x ^ — 0, but Im a 2 ^ 7^ 0, 
the misalignment concerns the other two principal axes 
perpendicular to the common axis 3. 
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TABLE I: The matrix of inertia Q: leading terms and relevant first-order terms. Other first-order terms are indicated with the 
symbol ~ 0. As the matrix is symmetric, first-order terms in the last three columns are not explicitly shown. Here J\ = 3(/3| + 
2/3 3 2 ) + 2^3(/3 2 2 72 + V5/? 3 2 73); Ji = 3(/3f+2/3f)- 2V3(/3| 7 2 +V5/Sg 73) ; and J 3 = A(/3hi+PH) + ^(X 2 +Y 2 ) + 2( V 2 + C 2 )+&e ■ 
The determinant of the matrix is G = DetC7 = 1152^f/3| + 2/3|) 2 (^f + 5/3|) _1 (/3| 7 2 + VE /3 2 7s) 2 • 
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C. The classical expression of the kinetic energy 

Now it is possible to express the classical kinetic energy 
(as given by Bohr hydrodynamical model) in terms of the 
new variables and of the intrinsic components q K of the 
angular velocity. The classical expression has the form 



(7) 



where Q = {£1, £ 2 , £9, qi, qi, 33}, ^ (m = 1, — , 9) are 
the time derivative of the nine parameters we have just 
defined, and q±, g 2 , 93 are the intrinsic components of the 
angular velocity of the intrinsic system with respect to 
an incrtial frame. The elements of the matrix Q (leading 
terms and relevant first-order terms) are shown in the 
Table H| 

The determinant G = Det Q takes the form 



Gcc/3 2 2 /3 3 2 (/? 2 2 + 2/?3 2 ) 2 (# 



5/? 3 2 ) 



L (/? 2 2 72 + V5/3 3 2 73) 



and, at the limit P3 <^ 02, results to be proportional to 
/3 2 ; and therefore consistent with that of the Bohr model 
for a pure quadrupole motion. This is a consequence of 
our choice of the normalization factors in the eq.s JSJ. 
This choice has other advantages: all the non diagonal 
terms involving the time derivatives of /3 2 or /3 3 and ei- 
ther the derivative of one of the other intrinsic amplitudes 
or </3 turn out to be zero in the present approximation. 
Other non diagonal elements are small (of the first or- 
der) in the "small" amplitudes 72, 73, X, Y, £, 77, (. In 
situations close to the axial symmetry, they have negligi- 
ble effect on the results (see Appendix IHl . with the only 
exception of elements of the last line and column. The 



latter, in fact, are still small of the first order, but must 
be compared with the diagonal clement ^3 , which is small 
of the second order in the "small" non-axial amplitudes. 
These terms play an important role in the treatment of 
the intrinsic component of the angular momentum along 
the approximate axial-symmetry axis, which will be dis- 
cussed in the next paragraph. 



D. Intrinsic components of the angular momentum 

According to the classical mechanics, the components 
L\, L2 and L 3 of the angular momentum in the intrinsic 
frame are obtained as the derivatives of the total kinetic 
energy with respect to the corresponding intrinsic com- 
ponent of the angular velocity: 

dT 

Lk = ^~. (8) 
oq k 

The part of the kinetic energy depending on the compo- 
nent qk has the form 

T k {qk) = \j k ql + F k q k (9) 

where F k is a function of the dynamical variables £j and 
of their derivatives with respect to the time, and is small 
of the first order according to our definition. As for the 
moments of inertia, 3% is small of the second order, while 
Ji, Ji are not small. According to eq.©, we have 



Lk — JkQk + Fk 
L k -Fk 



qk 



3k 



(10) 



4 



T k (L k ) 



T 2 



Yi (L t -F t ) 



Ik ~ F "' 



the substitutions 



2J k 2J k 



(11) 



For k = 1 , 2 the second term is small of the second order 
and can be neglected. It is not so for k = 3, as Jz is also 
small and of the same order as F%. In more detail, we 
have 



6(YX-XY) +2(Cry-7 ? C) 



where we have put 



7 



= V5 



72 - 73 



(12) 



(13) 



At this point, it will be convenient to express the vari- 
ables 72 and 73 as linear combinations of two new vari- 
ables, one of which is, obviously, 7 = \/572 — 73. The 
other one, that we call 70 , can be chosen proportional to 
the linear combination which enters in the expression of 
the determinant G, 



7o = c f/3f 72 + V5 ^373 
With this choice, we obtain 



(14) 



X = 

Y = 

V = 

C = 

£ = 

7 = 



wsint? 
wcosi) 
v simp 
v cos if 
wsinx 



V2 (^0. 



% + 50 2 / /3 2 /3 3 ] wcosx 



j = f(02,03)^0 2 2 +50 2 «o 



(17a) 
(17b) 
(17c) 
(17d) 
(17e) 

(17f) 
(17g) 



which gives for the determinant of the matrix Q 



ul v 2 u 2 w 2 f 2 (l3 2l f3 3 ). 



G = Dct Q = 2304 [0\ + 20. 

(18) 

The choice of the function f{02,0z) is irrelevant for 
what concerns the angular momentum. Non diagonal 
terms (small of the first order) would depend on this 
choice, but their effect is negligible (see Appendix[5J|. As 
a criterion to define the form of the function /, we observe 
that for permanent quadrupole deformation 02 — 02 and 
at the limit /3| <C 1 our value of G must agree with the 
result given, at this limit, by Eisenberg and Greiner [34| . 
This happens if the function / we have left undetermined 
tends to a constant when fa — ¥ 02- We adopt here the 
simplest possible choice, i.e. f [02,0s) = 1, to obtain the 
matrix Q given in the Table ITT1 and (at the leading order) 



72 = 



73 



7o/c- 



V5/? 3 2 7 



Pi 



5/31 



V5 70/c - (ih 



Pi 



5/31 



(15) 



and, at the leading order, 

PiPl 



Plil + Ptil 



Pi 



5Pi 



Pi 



5/3 3 2 



(7o/c) \ (16) 



In deriving the eq. 1161 the factor 1/c has been considered 
constant. We may note, however, that the same result 
holds at the leading order also if 1/c is a function of /3 2 
and/or fa- In fact, terms involving the time derivative of 
1/c also contain the "small" quantity 70, and their effect 
is negligible in the present approximation (see Appendix 
lB)l . E.g., one could choose for 1/c a quadratic expres- 
sion in P2, 03, in order to obtain for 70 an adimcnsional 
quantity, like 72, 73 and 7. The same argument applies 
for possible redefinitions of other "small" variables, like 



7- 



The expression (|12|) can be somewhat simplified with 



TABLE II: The matrix of inertia Q after the introduction of 
the new variables wo, v, u, w, tp, \, an d $ (see text). Here, 
Ji=j2= 3(/3| + 2j9f), and J 3 = 4 u 2 + 2v 2 + 8u 2 + 18w 2 . 
Only the leading terms are shown. Neglected terms are small 
of the first order (or smaller) in the sub-matrix involving only 
02, 03, Uo, v, u, ib, qi and q2\ of the third order (or smaller) 
in the sub-matrix involving only cp, x, $ and qs; of the second 
order (or smaller) in the rest of the matrix. 
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Det g = 2304(/? 2 + 2/? 2 ) 2 u 2 v 2 u 2 w 2 (19) 

(20) 



4ul + 2v 2 + 8u 2 + 18w 2 
Pi + Pi + ul 



+ 2(v 2 + v 2 tp 2 ) + 2(ii 2 + u 2 x 2 ) + 2(w 2 + w 2 ^ 2 ) 



2q 3 



2v 2 ip + Au 2 x + 6w ■& 



(21) 



E. 



Classification of elementary excitations with 
respect to K* 



We can now deduce the intrinsic components of the 
angular momentum: 



Li 
L 2 



Jiqi 

Jiq-i 
L 3 = J 3 q 3 + 



(22a) 
(22b) 



2v 2 ip + Au 2 x + Qw 2 -d 
2v 2 (q 3 + cp)+ Au 2 (2q 3 + X ) + 6u> 2 (3q 3 + #) 



+^u 2 q 3 



(22c) 



In the same way, we can obtain the classical moments 
conjugate to Xi $ an d f ( we observe that none of these 
variables appears in the expressions of G or J3): 
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2v 2 {tp + q 3 ) 
2u 2 (X + 2q 3 ) 
2w 2 (d + 3q 3 ) 



(22d) 
(22e) 
(22f) 



Now we can solve the system of equations (|22|) with re- 
spect to the variables q±, q 2 , q 3 , 'f , X, We obtain 



qi = LxUx (23a) 

q 2 = L 2 /J 2 (23b) 

q 3 = — (L 3 -p ip - 2p x - 3p$) (23c) 
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The equations (|23|l have a very simple meaning in the 
case where the potential energy does not depend on the 
variables ip, x or I n such a case (a sort of model ip~x~ 
^-instable, in the sense of the 7-instable model by Wilcts 
and Jean |3||) the conjugate moments of these three an- 
gular variables are constants of the motion, with integer 
eigenvalues n v , n x and n# (in units of Ti). Moreover, if 
we assume that uq — > 0, the third component q 3 of the 



angular velocity tends to 00 unless L 3 — p v — 2p x — ip^ = 
(cq. 123cj l). In this case, the operator L 3 is diagonal, with 
eigenvalues K = n v + 2n x + 3n$, and the three degrees 
of freedom corresponding to ip, x an( i ^ can be associ- 
ated to non-axial excitation modes with K =1,2 and 3, 
respectively. 

To investigate the character of the degree of freedom 
described by the parameter uq, and for a deeper under- 
standing of the nature of the other degrees of freedom, it 
is necessary to express the complete Hamiltonian in the 
frame of a definite model which, although not unique, is 
at least completely self-consistent at the limit close to the 
axial symmetry. 

In fact, it is now possible to use the Pauli prescriptions 
(37) to construct the quantum operator T corresponding 
to the classical kinetic energy T of eq. (|21l) . In doing 
this, we make use of the partial inversion of the matrix Q 
given by the solution (|23|l of the linear system (|22|l . and 
note that none of the variables ip, x, $ or 93 enter in the 
expression of G: 
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(24) 



This is, admittedly, only a semi-classical discussion. 
However, the formal quantum treatment in the frame 
of the Pauli procedure, shown in the Appendix IU1 gives 
exactly the same result. 

Until now, no assumption has been made on the form of 
the potential-energy operator, which will determine the 
particular model. A few general remarks on this subject 
are contained in the Appendix [D] We now assume that 
the potential energy can be separated in the sum of a 
term depending only on uq and another containing the 
other dynamical variables. In this case the differential 
equation in the variable uo is approximately decoupled 
from the rest. One obtains the Schrodinger equation 



1 d 

u du 
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L 2 [-^"o 
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U(u )} 



<j>(u o ) = (25) 



6 



where we have put 

Q, UQ = K — n v — 2n x — 3n§ . (26) 

If we assume, for simplicity, a harmonic form for the po- 
tential U(uo) — \Cu^, the eq. (|25l) is the radial equation 
of a bidimensional harmonic oscillator. For the existence 
of a solution, it is required that 

E U0 = {N U0 + l)huj U0 (27) 

with n Uo positive or negative integer and the integer 
N Uo > \n Uo \. Excitations in the degree of freedom cor- 
responding to the variable i*o carry, therefore, two units 
of angular momentum in the direction of the 3rd axis of 
the intrinsic reference frame. 

We could extend our model to include all the intrin- 
sic variables different from fa, fa- We assume a po- 
tential energy corresponding to the sum of independent 
harmonic potentials in the variables v, u, w and uq, plus 
a term depending on fa and fa (at the moment, we do 
not need to define the form of this term). We also as- 
sume that the equations in fa , fa can be approximately 
decoupled from those of the other variables and that, in 
the latter, fa and fa can be replaced by their average 
values. It is easy to verify that the differential equations 
in the pair of variables v, <p (or u, \ or corre- 
spond again to a bidimensional harmonic oscillator and 
that, as long as we neglect the rotation- vibration cou- 
pling, the eigenvalue K of the intrinsic component L3 of 
the angular momentum is given by 

K = ritp + 2n x + in^ + 2n Uo . (28) 

The energy eigenvalues are, for the equation in the vari- 
ables v, (p, 

E v = {N v + 1) huj v with N v > \n v \ (29) 

and have a similar form for the other two oscillators. 

It remains to consider the character of the different 
dynamical variables with respect to the parity operator. 
We know that the parity of the amplitude a^, is (— 1) A . 
Therefore, fa and 72 are even, while fa, X and Y are odd. 
As for 73, we observe that fa sin 73 is odd, and therefore 
73 must be even. As a consequence, are also even the lin- 
ear combinations 7 and 70/c defined in the eq. (|13I14[1 . 
The new variables i], £ and £ are odd, as (e.g.) fat; is an 
octupole amplitude and therefore is odd, while faS, must 
be even. Finally, on the basis of eq.s (|17|l we realize that 
v, u and w must be odd (while </?, x an d 1? are even). 
We have therefore identified elementary excitations cor- 
responding to K n = 1~ , 2~, 3~ and 2 + . Excitations with 
K = (of positive or negative parity) conserve the axial 
symmetry, and are related to the variables fa and fa. 
A particular example will be discussed in the following 
Sections. 



III. A SPECIFIC MODEL: AXIAL OCTUPOLE 
VIBRATIONS IN NUCLEI WITH PERMANENT 
QUADRUPOLE DEFORMATION 

Specific assumptions on the form of the potential- 
energy terms for all the variables describing the 
quadrupole and octupole degrees of freedom are neces- 
sary in order to obtain definite predictions, also if these 
are limited to the axial modes. 

We discuss here, as an example, the case of axial oc- 
tupole excitations in nuclei which already possess a stable 
quadrupole deformation. In this case, one obtains rela- 
tively simple results, suitable for comparison with exper- 
imental data. This comparison will be performed in the 
next Section. 

Following the usual treatment of vibration + rotation, 
we put fa = fa + P' 2 , with fa = constant and \fa 2 \ ^ 
The new variable j3' 2 is therefore assumed to be small (of 
the first order) as all other variables, with the exception 
of fa. With this choice, and assuming that the variables 
introduced in the eq. (f 1 71) are suitable to describe the 
other degrees of freedom, the matrix Q takes a form sim- 
ilar to that of Table II (with fa substituted by fa and fa 
by $ 2 ) and ~ a t the lowest significant order - results to be 
diagonal with respect to the variables fa, /3 2 , <7i and q 2 . 
Moreover, in our model, the amplitude of oscillation for 
all degrees of freedom different from fa are constrained 
to very small values: this fact implies strong restoring 
forces and, therefore, oscillation frequencies much larger 
than for ^3. 

In the limit of small amplitude of the octupole oscil- 
lations, this case has been discussed, e.g., by Eisenberg 
and Greiner [34|. In their approach, the "intrinsic" ref- 
erence frame is chosen to coincide with the principal axis 
of the quadrupole deformation tensor, but the differences 
between their approach and ours tend to disappear for 
\fa\ <C fa- The model we try to develop should be able 
to describe (axial) octupole vibrations of finite ampli- 
tude, but its limit for \fa\ <§; fa must obviously agree 
with the results of Eisenberg and Greiner. 

The quantum-mechanical equation of motion for fa 
can be obtained with the Pauli prescription 37], with 
the additional assumption that the equations involving 
fa and P' 2 or the angular-momentum components L±, L 2 
are effectively decoupled from those containing the other 
dynamical variables and/or the L3 operator. The latter 
equations could possibly be complicated, and substan- 
tially coupled with one another and with the angular 
momentum component L3 along the (approximate) sym- 
metry axis. A short discussion of this subject, with some 
simplifying assumptions, has been given in the previous 
section Til Dl At the moment, we assume that terms in- 
volving 02, Li and other dynamical variables different 
from fa contribute to the total energy with their own 
eigenvalue, independent from the eigenfunction in the fa 
degree of freedom, and we only consider their lowest- 
energy state. We also assume that this state has K = 
(and neglect, as usual, the possible rotation- vibration 
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coupling) 
the form 



In this case the complete wavefunction has 



* oc ^(^ 3 )$o(^,7o,^, i o,eo ) ^X I <PpM,o 

cx iK/? 3 )$o Yj M (£l) (30) 

where V J M M , are the Wigner matrices and (3 3 ) — 
(—) J ip(f3 3 ). The differential equation for the wavefunc- 
tion ip{(3 3 ), obtained with the Pauli prescription, has the 
form 



/--1/2 d 1/2 d 
d/3 3 d/3 3 

h 2 J(J + l) 



+ vw 

v(/3 3 ) = (3i) 



where V(/?3) is the potential-energy term. The expres- 
sion of the determinant G is not uniquely defined, as 
it also depends on the part of the matrix of inertia Q 
involving all other dynamical variables. This is a gen- 
eral problem for all models where part of the degrees of 
freedom are ignored (and, in collective models, a num- 
ber of dynamical variables describing the details of nu- 
cleon degrees of freedom, are certainly ignored). With 
the present choice of dynamical variables (Table [HJ, one 
gets G =Det Q = 2304 u% v 2 u 2 w 2 (0 2 + 2f3f) 2 , and the 
eq. <j3*Tjl becomes 



d 2 



4/3 3 



+ 



h 2 j(j + l) 

2 Ji 



2/? 3 2 d(3 3 



+ Vtfa) (32) 



V>(/3 3 ) = £Ws) ■ (33) 



This is the equation we have used in ref.s 0,^3- A differ- 
ent choice of the dynamical variables would have brought 
to a different result. E.g., with the variables used in the 
Table [I] one obtains 

G cx p 2 0lC0 2 2 + 20l) 2 CPi + 5(3 2 )-Hp 2 l2 +_V5(3 2 l3 ) 2 . 
However, in this case the limit \f3 3 \ <C |/?2| would not 
correspond to the Eisenberg-Greiner result, due to the 
presence of the factor /3f in the expression of the de- 
terminant G. In fact, when the spherical symmetry is 
broken by the permanent quadrupole deformation, and 
at the limit of small octupole deformation, the octupole 

amplitudes aji are decoupled from one another [34| and 
it would have been more reasonable to chose a dynam- 
ical variable u 3 = ^373 in the place of 73. With this 
substitution, the factor /3| disappears. 

It is convenient to express the differential equation 133|) 
in terms of the adimensional quantities x = V2/33//?2 and 
e = j^PlE, v = j^fl^V. One obtains 



d 2 ip(x) 2x d^(x) 



dx 2 



1 + x 2 dx 



JQ7 + 1) 
6(1 + a; 2 ) 



v(x) 



ip(x) = 
(34) 

where ip(—x) — (—) J ip(x). This equation reduces to that 
of Eisenberg and Greiner [34| when x <C 1. 



As for the potential v(x), we have explored two possible 
forms: a quadratic term 1 v(x) — ^cx 2 or a square-well 
potential, as it has been adopted [2] at the critical point 
in the X(5) model ( v(x) — for \x\ < b and = +00 
outside, so that ij)(±b) = ). In both cases, there is 
a free parameter to be determined from the comparison 
with empirical data. 

We now discuss in particular the second case. We may 
note that for V = the eq. I|34|l is formally equivalent to 
that of spheroidal oblate wavefunctions (see eq. 21.6.3 
of ref. 38]) with the parameters to, A and C redefined as 
to = 0, C 2 = e and A = J (J + l)/6 — e. Here, however, 
the solution is confined in the interval —b < x < b and 
the equation has been solved numerically. For a given 




FIG. 1: Examples of wavefunctions ip 3 j(x) as a function 
of x/b, for b — 1.73. Pairs of curves refer to consecutive 
values of J, with even parity (continuous line) or odd parity 
(dotted line). Part a : ,r = 0+,l", b : J n = 10+, 11", 
c : J n = 18+, 19" with s = 1; d : J™ = + , 1" with s = 2. 
With increasing angular momentum, the difference in |i/)(a;)| 2 
between consecutive values of J tends to vanish and, as a 
consequence, the positive and negative part of the band merge 
together. 

b and for every value of J, one obtains a complete set 
of orthogonal eigenfunctions, with an integrating factor 
(1 + x 2 ). These eigenfunctions can be characterized by 
the quantum number s = v + 1, where v is the number of 
zeroes in the open interval < x < b. A few examples of 
wavefunctions corresponding to the square-well potential 
with b — 1.73 and for different values of s and J are 
depicted in the Fig. ^ 

The dependence of the eigenvalues on the parameter b 
is illustrated by the Fig. |3 where the ratios E(J)/E(2) 
are shown for the g.s. band (s = 1). Other possible 
choices of the set of independent dynamical variables 
would have brought to a different equation, but the dif- 
ference would have concerned the coefficient of the first- 



1 After completion of this work, we have been informed that a 
quadratic potential plus a centrifugal term with variable moment 
of i ner tia has also been considered in the model by Minkov et 
al. El. 
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FIG. 2: (Color online) Ratio E(J n )/E(2 + ) as a function of 
b, for states of the ground-state (s = 1) band with different 
J". Dotted lines show the results obtained with a differential 
equation corresponding to that of Eq. 1341 but without the 
Brst-derivative term. 



derivative term the eq. (|34[1 , with very small effect on the 
results, as long as b is in a range of "reasonable" values. 
To exemplify the effect of this term, results obtained with 
the coefficient of the first derivative put to zero are also 
shown - as dotted lines - in the Figure |21 Differences 
between the two sets of results turn out to be very lim- 
ited for small values of b (at least, up to b w 2). There 
is at least one case ( 226 Th) in which our results with 
the square-well potential are in good agreement with the 
level scheme, for low-lying states of positive and negative 
parity, while for 228 Th a better agreement is obtained 
with the quadratic potential. The possible interpreta- 
tion of this result as evidence for a phase transition in 
the octupole degree of freedom is discussed in the next 
Section. 



IV. EVIDENCE OF PHASE TRANSITION IN 
THE OCTUPOLE DEGREE OF FREEDOM 

A. The Radium and Thorium isotopic chain 

A phase transition in the nuclear shape manifests itself 
as a relatively sharp change of a proper order parameter 
- e.g., the ratio R = E(4 + )/E(2 + ) - as a function of a 
driving parameter which can be, in our case, the number 
of neutrons in the isotopes of a given element or the num- 
ber of protons along an isotone chain. Due to the finite 
number of degrees of freedom, the transition region has 
a finite width around the critical point, and extends over 
several nuclides in the chain. In the case of transitions be- 
tween spherical shape and axial quadrupole deformation, 
the X(5) symmetry, valid at the critical point, predicts 
E(4 + )/E(2 + ) = 2.91 and we can use this criterion to 
locate the critical point of the phase transition. The sit- 
uation is more complex when the quadrupole and the oc- 
tupole degrees of freedom must be considered at the same 
time. Fig.|3|shows, as a function of the neutron number of 
Ra and Th isotopes, a few parameters which can be used 




- 220 230 A Th 

\K > 



A A 



if = 



d 



130 140 N 

FIG. 3: (Color online) Indicators of the quadrupole collectiv- 
ity (left) and of the octupole collectivity (right), as a function 
of the neutron number N in the isotopic chain of Ra (circles) 
and Th (triangles) : a - Excitation energy of the first 2 + level; 
b - Energy ratio E(A + ) / E(2 + )\ c - Excitation energy of the 
first level of the K n — CP band, Jq = 1~ (open symbols) 
and of the lowest known level of other negative-parity bands, 
J f 7 = 2" or Ij (full symbols); d - Energy ratio E(\-)/E{2+). 
The horizontal line in the part b shows the value (2.91) ex- 
pected for the X(5) symmetry. Data are from ref.s |39ll40|. 



as indicators of the quadrupole and octupole collectivity. 
As for the quadrupole mode, the decrease of E(2 + ) with 
increasing N (Fig. |3Ji) shows a corresponding increase 
of collectivity. Moreover, in the Fig. [3f> we observe the 
transition between the vibrational (or not collective) be- 
havior of the lighter isotopes of the chain and a clear rota- 
tional behavior (R w 10/3) above A = 226, with a critical 
point which can be located around A = 224. The ratio 
E{1~)/ E(2 + ), depicted in Fig.|3Jf, shows that the relative 
importance of the octupole collectivity increases with de- 
creasing TV and reaches its maximum in the region below 
N = 138, where the critical point of the phase transition 
in the quadrupole mode could be located on the basis of 
Fig. [3J>. Heavier isotopes show evidence of octupole vi- 
brations (of different K) around a quadrupole-deformed 
core 0,H2|. Lighter isotopes (N < 132) appear not to 
be deformed in their lower-J states. However, at larger 
angular momentum a rotational- like band develops, and 
this band has the alternate-parity pattern typical of a 
stable octupole deformation |43| . 

The model introduced in the first part of this paper, 
and developed in the section ITTT1 for the particular case of 
a permanent quadrupole deformation, assumes that non- 
axial amplitudes are constrained to very small values by 
the large restoring forces. This implies that excitation 
of one of the non-axial degrees of freedom leads to high 
excitation energy, compared to that of the first 1~ level of 
the K = band. Experimental data of Fig. show that 
this is actually the case for the light Thorium isotopes, at 
least up to A — 228. In fact, the first 1~ level is not far 
from the first 2 + and much lower than levels belonging 
to negative-parity bands with K ^ (as the lowest 2~ 
or the second 
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FIG. 4: (Color online) Ratios E(J W )/E(2 + ) as a function of J, for positive-parity states (circles) and negative-parity states 
(triangles) of the ground-state (s = 1) band of 226 Th and 228 Th, compared with different model calculations: rigid rotor (curve 
a, even parity only), present model at the critical point (curves b and b'), present model with harmonic potential in (curve 
c). The curves b and c correspond to a fit on the lowest V state, the curve b' , on the 20 + state. Note that the even-parity 
parts of the curves b and b' are very close to the curve a for J < 6, while, for each one of the curves b, b' and c, the even- and 
odd-parity branches tend to merge together at large values of J. 



B. Comparison with experimental data for ' Th 

Our model assumes a permanent quadrupole deforma- 
tion. Therefore, it can be useful only for relatively heavy 
Th isotopes (Fig. Eh, b). The quadrupole-deformed re- 
gion extends above the mass 224 (which could corre- 
spond to the critical point of the phase transition, hav- 
ing E(A + )/E{2 + ) fn 2.91). Heavier Th isotopes (with 
A > 230) show negative-parity bands built on the differ- 
ent states of octupole vibration, from K = to K = 3, 
with band heads much higher than the first 2 + (Fig. • 
Only for lower A, the 1~ band head of the K v — 0~ band 
decreases well below the band heads of all other octupole 
bands, and higher levels of the K* — CP band merge with 
those of positive-parity of the ground-state band (Fig.^J, 
approaching (but not reaching) the pattern expected for 
a rigid, reflection asymmetric rotor. 

The region of possible validity of our model is therefore 
restricted to 226 Th and 228 Th. 

The ratios E{J 7T )/E{2+) for the low-lying states of 
226 Th and 228 Th are depicted in the Fig. g| and com- 
pared with the predictions of different models. The rigid- 
rotor model cannot account for the position of the lowest 
negative-parity levels, and overestimates the excitation 
energy for all the high-spin states. Instead, a rather good 
agreement is obtained with the present model, if one as- 
sumes, in the case of 226 Th, a square-well potential (as 
the one hypothesized by Iachcllo in his X(5) model) and, 
in the case of 228 Th, a harmonic restoring force. In both 
cases, the free parameter of the model has been adjusted 
to reproduce the position of the 1~ level. As shown in the 
Fig. 21 (dotted line) a much better agreement with the 
high-spin levels of 226 Th is obtained with a slightly dif- 
ferent value of the parameter (b = 1.87 instead of 1.73), 



at the expense of a very limited discrepancy for the l - 
level. 

Therefore, for what concerns the level energies of 
the ground-state band (including in it also the odd- J, 
negative-parity states), 226 Th seems to present the ex- 
pected behavior of a nucleus with permanent quadrupole 
deformation and close to the critical point of the phase 
transition in the octupole mode. 



C. Other possible tests of the critical-point 
behavior 

A considerable amount of experimental information 
has been reported, in the last few years, on possible can- 
didates 0, IE 0, i, H [H IH E Hi for the dynamical 
symmetry X(5) (phase transition point in the quadrupole 
mode). It is now clear that the agreement between ex- 
perimental and calculated energies for the ground-state 
band does not automatically imply that a similar agree- 
ment exists also for other observables, like the excitation 
energy of the second + level (the band head of the s = 2 
band) and the in-band and inter-band transition proba- 
bilities. In several transitional nuclei, the excitation ener- 
gies in the ground-state band are in excellent agreement 
with the X(5) predictions but the calculated ratios of the 
B(E2) transition probabilities fail to reproduce the ex- 
perimental ones 0, unless an ad-hoc second-order term 
is included in the E2 transition operator [l4L |44| . It is 
therefore important to test the predictions of our models 
also for what concerns such observables. 

The low-lying level scheme of 226 Th is shown in the 
Fig. [5J together with the one resulting from the present 
model, with the value of b adjusted in order to reproduce 
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FIG. 5: Experimental level scheme of Th compared with 
the predictions of the present model at the critical point. The 
two spectra are normalized on the first 2+ level, and the 
model parameter b is adjusted to reproduce the position of 
the first 1". 



TABLE IV: Experimental and calculated values of the ra- 
tios of reduced strengths (in W.u.) R = B(E1)/B(E2), for 
transitions coming from the same level of 226 Th. 



Trans. 1 Trans. 2 



R x 10 5 



Ji 






Theoretical 


Experimental 


8+ 


El 


7~ 


E2 


6+ 


1.3 


2.0 


(8) 


9" 


El 


8+ 


E2 


7" 


1.3 


1.7 


(2) 


10+ 


El 


9~ 


E2 


0+ 


1.6 


1.5 


(I) 


11" 


El 


10+ 


E2 


9" 


1.6 


1.7 




12+ 


El 


ll - 


E2 


10+ 


1.8 


1.6 


(i) 


13" 


El 


12+ 


E2 


11~ 


1.8 






14+ 


El 


13- 


E2 


12+ 


1.9 


1.4 


(i) 


15" 


El 


14+ 


E2 


13" 


2.0 


1.7 


(3) 


16+ 


El 


15" 


E2 


14+ 


2.1 






17" 


El 


16+ 


E2 


15" 


2.1 


1.5 


(3) 


18+ 


El 


17" 


E2 


16+ 


2.2 






19" 


El 


18+ 


E2 


17- 


2.3 


1.7 


(4) 


"Values used for normalization. 



the empirical value of the ratio E(l~)/E(2 + ). At the mo- 
ment, only the first two levels of the s — 2 band (0^, 2^~) 
are known and their excitation energies are somewhat 
higher than the values predicted at the critical point. 
We can observe, however, that also in the best X(5) nu- 
clei 0, 03, the position of the levels of the s = 2 band 
deviates somewhat from the model predictions (although 
in the opposite direction). In our opinion, a similar qual- 
itative agreement is obtained also in the present case. 
The negative-parity levels of the s = 2 band are pre- 
dicted to lie at higher energies, and could be difficult to 
observe. Absolute values of the transition strengths are 
not available for 226 Th, but some relevant information 
is provided by the branching ratios in the level decays. 
In the tables IIIIIIVI experimental ratios of the reduced 
transition strengths for El or E2 transitions coming from 
the same level are compared with the model prediction 
at the critical point. 



TABLE III: Experimental and calculated values of the ratios 
of reduced strengths for El transitions coming from the same 
level of 226 Th. 





Trans. 1 


Trans. 2 


B^ h (E1)/S^ /2 (E1) 








Theoretical Experimental 


i" 


^0+ 


^2+ 


0.47 0.54 (5) 


3- 


^2+ 


^4+ 


0.65 0.99 (25) 


2+ 


=> 1" 


=> 3" 


0.63 0.60 (18) 



The electric dipole moment would vanish for the col- 
lective motion of a fluid with uniform charge density, 
as the center of charge would coincide with the center 
of mass. Therefore, the observed El transition ampli- 



tudes are entirely due to the non uniformity of the nu- 
clear charge distribution [45|. To calculate the value of 
B(E1)= (i\\M(El)\\f) 2 /(2J l + 1), the El transition op- 
erator has been assumed to have the form |46l I47L l48l |49| 

M^El) = CifolhYi,,, , (35) 

with the constant factor C\ depending on the nuclear 
charge polarizability. The E2 transition operator for the 
in-band transition and at the limit close to the axial sym- 
metry has been taken in the simple form 

A/( M (£2) = C 2 (/3 )Y 2 ^ , (36) 

neglecting the (weak) dependence on Therefore, the 
theoretical ratio of the reduced strengths for transitions 
of different multipolarity (El and E2) is determined apart 
from a constant factor, which must be fixed by compar- 
ison with the experimental data. The average of the ra- 
tios J5(E1)/B(E2) for the transitions coming from the 
10+ and 11" levels has been used for normalization of 
the theoretical values given in the Table IIVI 



V. CONCLUSIONS 

A theoretical scheme for the description of quadrupole 
plus octupole excitations close to the axial symmetry 
limit has been developed in the Section [H] and special- 
ized in the Section 1 1 1 1 1 to the simpler case of a perma- 
nent (and axially symmetric) quadrupole deformation. 
In principle, the model should be able to describe the 
wide field of reflection-asymmetric nuclear shapes, close 
to (but not coincident with!) the axial symmetry limit. 
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Calculations of nuclear shapes in the frame of the HBF- 
Cranking model |50| for nuclei of the Radium - Thorium 
region find a large variety of results, including proper 
potentials for quadrupole-octupole vibrations around a 
spherical shape or for octupole vibrations around a de- 
formed, reflection symmetric shape, and also situations 
with a rather flat minimum of the potential along a line 
at constant f3 2 with \0a | < (3% iax . The latter case is just 
what is expected for our "critical point" of the phase 
transition in the octupole mode. In the Section II VI we 
have investigated the evolution of the nuclear shape along 
the isotopic chain of Thorium, and shown that evidence 
of phase transition exists, not only in the quadrupole 
mode but also in the octupole mode around a stable 
quadrupole deformation. The model developed in the 
Section ITTT1 results to be able to account for the experi- 
mental data of 226 Th (at the critical point of the phase 
transition in the octupole mode), and also of its neighbor 
228 Th (characterized by axial octupole vibrations). More 
and improved experimental data on E2 and El transition 
strengths would be necessary for a more stringent test of 
the model predictions. 

Further developments of the calculations are in pro- 
gram, to provide detailed predictions for other signifi- 
cant cases: first of all, those of 224 Ra and 224 Th, whose 
positive-parity levels show an energy sequence very close 
to Iachello X(5) predictions for the critical point of phase 
transition between spherical shape and axially symmetric 
quadrupole deformation. 



APPENDIX A: SUMMARY OF THE GENERAL 
FORMALISM 

The general formalism to describe collective states of 
rotation/ vibration in nuclei is discussed e.g. in ref. |18| . 
The nuclear surface is described in polar coordinates as 



1+ E E < x) nA e ^) 

A /j=-A,A 



(Al) 



In the sum, the values of A are now limited to 2,3. A 
term with A = 1 should be included in order to maintain 
fixed the position of the center of mass |34| : 




^v/(2A + l)(2A' + l) 



Q 



(A) 



a 



(A') 



(1) 
I" 



(A2) 



This term, however, is inessential for the following dis- 
cussion, and will be omitted. For not-too-large deforma- 
tion, in fact, the amplitudes otjP are much smaller than 
the others, and their effect results to be negligible for 
most purposes (with the noticeable exception of the El 



transition amplitudes). At the same level of accuracy, 
one can also neglect the slight variation of Ro necessary 
to keep the volume exactly constant. 

In the Bohr model that we are considering here, the 
classical expression of the collective kinetic energy is 



T 



,(2) 



2 3 



E 

^— — 3,3 



a 



(3) 



(A3) 



In order to simplify the notation in the following, we 
include the inertia coefficient B\ in the definition of the 
collective variables a]^ . In the literature, this symbol is 
usually reserved to the variables defined in the intrinsic 
reference frame and, from now on, we will always use this 



reference frame for the collective variables a 



B\ a 



(A) - 



E< 



.(A) jD (A)* 



(A4) 



where are the Wigner matrices and 9i the Euler an- 
gles. The axis of the intrinsic reference frame are defined 
along the principal axis of the inertia tensor. 

The expression of the kinetic energy in terms of the 
time derivatives of the intrinsic deformation variables 
and of the angular velocity q of the intrinsic frame 
with resp ect to an inertial frame is discussed, e.g., in the 
ref.s |l8tl34l| . If only the quadrupole mode is considered, 
the total kinetic energy can be expressed as the sum of 
a vibrational term (in the intrinsic frame) and a rota- 
tional term. If the octupole mode is also considered, a 
rotation- vibration coupling term must be added |34l l35| : 



T = T, 



vib "i~ Trot 



coup 



where 



(A) 

A, 



(A5) 
(A6) 



= - ^ Qkfjk'Jkk' 



k,k' 



(A7) 



k,k' 



E * 



E\/ 3 ( 2A+1 ) 



K } ) , 



(i) 



(0) 



(A8) 



Here, % (k = 1, 2, 3) are the Cartesian components of the 
angular velocity q along the axes of the intrinsic frame, 
while the are (2A + l)-dimensional matrices giving 

the quantum-mechanical representation of the Cartesian 
components of an angular momentum M = A in the in- 
trinsic frame, subject to commutation rules of the form 

M\ X} M^ } - A'4 x h4 x} = -iM^ x} . (A9) 
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We assume, as usually, 
( M^) vu , = v5 v , v . 
( ) vv , = I (V(A-«/) (A+V+l) 

+ \/(A + j/) (A - 1/+ 1) <5„>_i) 

( M 2 {A} - ^ V (A + ^ + l) <W +1 
+ i V(A + ^) (A-I/ + 1) Ou> : i>— 1 I . 

Taking into account the properties of the M^ X \ it is pos- 
sible to obtain the explicit expression for the diagonal and 
non diagonal elements of the tensor of inertia [18L l35| 

J* = \ E «W (m^m^ + m^m^)^ . 

(A10) 

The spin operators transform as the Cartesian 

components of a vector under rotation in the or- 

dinary space. Taking into account the commutation rule 
(IA9|) , we can define the irreducible tensor components of 

as 



and that 



M, 



(1){A} 




M. 



(All) 



and express the products of two Cartesian components as 
the sum of products of two tensor components. To this 
purpose, we define, for each value of A, the irreducible 
tensor product 



T {j) 



(J) 



\m {1) <g> M {1) 

L J m 

Y^{^\ V '\Jm)M^ M ( J ] (A12) 



(where the common suffix {A} has been dropped, for sake 
of simplicity). We now introduce the reduced matrix 
elements of the tensor operator T^ J \ 

(A||tW||A) = (-1)V2J + t{* 1 4a(A+1)(2A + 1) 



to obtain 



Jrn 

= E 

Jrn 

X ; 

V2J+ 1 
We now observe that 



(A13) 

E(lM^V^) (tL J) ) , (A14) 

Jrn 

^(l M l M '|Jm)(A||TW||A) 

(A^A — v | Jro) . 



(A15) 



(J) 



J] a^a^ (Ai/A - i/| Jm) (A16) 



to obtain 



){A} 



J2^=^=(l^'\Jm) (A||T«||A) 



J,m 

x [ a ( A W A > 



{J) 



(A17) 



The rank of the tensor product of the two identical 
vectors M' 1 "*' must be even, and therefore the possible 
values of J are limited to and 2. Now we can substitute, 
in the eq. (|A10J) . the cartesian components of the angular 
momentum with its tensor components defined in the eq.s 
(|A11(1 . The possible values of m (and J) contributing to 
the sum are limited to m = (and therefore J =0 or 2) 
in the case of (M 3 ) 2 , to m = or 2 (J =0 or 2) for (Mi) 2 , 
(M 2 ) 2 and (MiM 2 ), to m = 1 (J = 2) for (M 3 Mi) and 
(M 3 M 2 ). One obtains 



ji = E i c °w 

A 
1 



a (A) a 



(A) 



(0) 




C 2 (A)Re 



a< A > ® a< A > 



(2) 


(2) 
2 



(Al8a) 



+ 



J2 = E C "( A ) 

A 

1 =C a (A) 



a< A > ® a< A > 



(o) 





V6 

+ C 2 (A)Rc 



«W ® a< A ) 



a< A > ® a< A > 



(2) 


(2) 
2 



(A18b) 



a< A > ® a™ 



(o) 





I C 2 (A) 







a( A ) ® a( A )" 







J7l2 
i7l3 

1/23 



EC2(A) Im [a (A) ® a (A) 

A 

E C 2(A) Rc 



where 



Co (A) 
C 2 (A) 



= E^W Im ® ° (A) 

A 

H) A ^V2ATT 



(2) 
2 

(2) 
1 

(2) 
1 



(Al8c) 
(A18d) 

(Al8e) 

(A18f) 



(-1) 



A+l 



A(A+1)(2A+3)(4A 2 - 1) 



30 
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and therefore 

C (2) = 2 ^5 , 
C 2 (2) = -v^T ; 



Co (3) = -4 V7 
C 2 (3) = 



14 



(A19) 



If we chose as the intrinsic reference frame the principal 
axis of the tensor of inertia, we must put Jn = Jn = 
1/23 = 0. Using Eq.s (11121-ill we obtain, up to the first 
order 



J12 = -2 V21 Im I a (2) <g> 5 (2) 

+ 6 V14 Im 

Jl3 + %Jt& = -2 V21 

+ 6 VU 



(2) 
2 



(2) 
2 



a< 2 > ® a< 2 > 



(2) 



= 



0, (A20a) 

(2) 

1 

(A20b) 



Here we make use of the fact that the zero-order terms 
are automatically set to zero if the ajL are defined ac- 
cording to eq.s (11121) . By inserting these definitions in 
the above equations, and retaining only the first-order 
terms, one obtains 



Im 



Im 



S« ® a« 



® 5< 3 > 



® a< 3 > 



(2) 
1 

(2) 
2 

(2) 
1 

(2) 



(2,0,2,112,1) A a?> 

1 R 

7j^P2 a l , 



(2,0,2,2|2,2)/3 2 Ima 



j (3 2 Im a 2 2) , 



(3,0,3,l|2,l)/? 3 af 5 



s( 2 ) 
2 



(3,0,3,2|2,2) /3 3 Im a£ 
J A Im 4 3) 



We obtain, therefore, from eq.s l|A20|) 

/? 2 ai 2) = -\/2/?3^ 3) , 
/3 2 Im 4 2) = -V5 /3 3 Im a£°- 



(A21) 



APPENDIX B: EFFECT OF NON DIAGONAL 
(FIRST-ORDER) TERMS 

The matrix Q, as it results from the Table [I] contains 
zero-order terms only in its principal diagonal (whose 
last element, however, is small of the second order). Non 
diagonal terms have been expanded in series up to the 
first order in the "small" dynamical variables (all of them, 
apart from /3 2 and /3 3 ). We shall call a generic "small" 
variable, different from (3\ (A = 2, 3). 



Terms of the last line and column (those related to the 
third intrinsic component of the angular velocity) are dis- 
cussed in the Section Til Dl Here, we consider a simpler 
problem: the inversion of a matrix Q which has finite 
values for all terms in the principal diagonal, and only 
"small" values for all others. In the zero-order approx- 
imation, the inverse matrix A = Q is diagonal, with 
diagonal elements = 1/G W . 

The first-order approximation gives the non-diagonal 
elements 



4 — A — — 



G 

G ^pG vv 



(Bl) 



We are interested in particular on the effect of non di- 
agonal terms on the coefficients of the derivatives with 
respect to [3% or /3 3 . 

Terms involving the derivatives with respect to [3\ and 
£/c have the form 



G~ 



d 



d(3x 



G~ 



d 

df3x 



G*A 



\1: 



(B2) 



The non diagonal matrix elements of the matrix A are 
of the first order in the "small" variables. They have 
therefore the form f^.(02, where fxk(02,0a) = 
dA\,k/ d£,i, or - possibly - are the sum of several terms 
like that. 

Substituting this expression in the Eq. (|B2|I one ob- 
tains 



G 



d 



G? A 



kX 



d 
Wx 



G~ 



G~-A 



d 
Wx 

w) +5tk 



Xk 




d 

W 

d 
Wx 
8 



(B3) 



The last line of Eq. (|B3|) only contains the partial- 
derivative operator with respect to the small variable . 
Compared with the diagonal term involving the corre- 
sponding second-derivative operator, it contains a small 
factor & more, and can be neglected. 

In the first line of the expression, all terms are po- 
tentially of the same order of magnitude of the leading 
ones, and in principle could not be neglected. The first 
of the terms in the square brackets contains the partial 
derivative with respect to In the spirit of the adia- 
batic approach, we try to estimate the expectation value 
of this term for the ground-state wavefunction in all the 
variables, with the exception of 02 and /3 3 . If we assume 
that the potential is harmonic for the ensemble of these 
variables, the ground-state wavefunction can be expected 
to be close to a multivariate Gaussian function. If, more- 
over, £fc and £j are uncorrelated, the expectation value of 
iid/dik is zero, as (&) = 0. 
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For the case i = k, instead, the expectation value is 



the corresponding quantum operator has the form 



/ exp(-< 2 /2)^ exp(-a£ 2 /2)d£ 



J exp(— a£ 2 )d£ 
Jexp(-aC 2 /2) (™2a£ 2 ) exp(-af 2 /2)d£ 



-1 



/ exp(— a£ 2 )d£ 



(B4) 



and (in this approximation) cancels the third term, Sik- 
In conclusion, the coefficient of d/d(3\ coming from the 
non-diagonal terms of the matrix Q can be approximated 
with a sum of expressions like 



/S(A,A) (| f ) 



(B5) 



which approximately vanishes for i ^ k, and also van- 
ishes for i = k unless both A\k and G depend explicitly 
on £fc- Moreover, also in this latter case, it is possible to 
eliminate non-diagonal elements of the matrix Q of the 
form Qxk = 5Afc(/32, ^3) £fc, with a slight change in the 
definition of (3\, without any other effect at the present 
order of approximation. Namely, it is sufficient to sub- 
stitute the dynamical variable (3\ with the new variable 



f- -— G~ 1/2 V— G 1/2 A — 
2 G A ^d^ 



(C2) 



where G = Det Q and A = Q~ x . The choice of the 
"best set" of dynamical variables is, in part, related to 
the expression of the potential-energy term, and it is not 
obvious that it will eventually coincide with the one dis- 
cussed in the section III El However, it can be useful to 
explore the properties of the kinetic-energy operator in 
the particular model in which the matrix of coefficients 
Q is exactly that of the Table [HJ with the non-diagonal 
terms confined in one single line (and column) of the 
lower 6x6 sub-matrix. 

From now on, we limit our discussion to the lowest 
6x6 sub-matrix Q\ of the matrix Q (and of the matrices 
derived from Q, Tables IVT1 IVII|I as the six corresponding 
variables-^, %, 1?, 6\, 8 2 , 83 - are effectively decoupled 
from the others. 

The formal procedure is discussed, e.g., in the Ch.s 5 
and 6 of the ref. [34]] • The first pass is the substitution 
of the intrinsic components q\ , #2, 93 of the angular ve- 
locity with the time derivatives of the three Euler angles 



(B6) 



to obtain 

fa = fa ~ gxk(Jh,lh)tk£k + 0(&) (B7) 
and therefore, up to the second order in 
P 2 +2s Afc (ft, 

= (fa -SA*(#!,ft)&&) 2 + 29Xk(p2,fo)Zkl3x£ 

» (fa)\ (B8) 

This argument can be easily extended to the case in 
which the dependence of A\k on comes from the de- 
pendence on of one of the terms at the denominator 
in the eq. QBlfl. 



APPENDIX C: QUANTIZATION ACCORDING 
TO THE PAULI RULE 

In this Appendix, we discuss some aspects of the quan- 
tization of the kinetic energy expression given, e.g., in 
the eq. filjl . by means of the Pauli procedure: namely, if 
the classical expression of the kinetic energy in terms of 
the time derivatives of the dynamical variables £ M is 



(CI) 



(C3) 



with 



V = 



— cos 8 3 sin 82 sin 83 
cos 6*3 
1 



sin 83 sin 82 
cos 82 



(C4) 



As a consequence of this substitution, the 6x6 matrix 
Qi transforms according to the relation 



g[ = w & w, 



w 



{1} {0} 

{0} v 



(C5) 



(where W is the transpose of the matrix W and {1}, {0} 
are the 3x3 unit matrix and null matrix) and takes the 
form shown in the Table The next step is the 

inversion of the matrix Q[. In doing this, the explicit 
form (|2(Jfl of J% has been introduced and second-order 
terms in the small amplitudes have been neglected. At 
this point, it is necessary to introduce the intrinsic com- 
ponents of the angular momentum operator [Li, L2, L3) 
in the place of the derivatives with respect to the Euler 
angles. The expression of in terms of d/ddk, and vice- 
versa are given, e.g. in the Ch. 5 of ref |3J|. One gets 



—2 



d8 k 



(C6) 
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TABLE V: The relevant part of the matrix of inertia, Q[, after the substitution of the angular- velocity components q k 
(k = 1, 2, 3) with the time derivative of the Euler angles 6 k . Here, Ji = J 2 = 3(/3| + 2/3|), and J 3 = 4 u% + 2v 2 + 8u 2 + 18w 2 . 
The determinant of this 6x6 matrix is now G[ = Det G[ = 288 ul v 2 u 2 w 2 (/3 2 + 2/3 2 ) 2 sin 2 6 2 ■ 







X 


1? 


e\ 


0-2 


83 




2v 2 








2v 2 cos 2 





2v 2 


X 





2u 2 





4u 2 cos 02 





4u 2 


d 








2w 2 


6w 2 COS 02 





Gw 2 


&i 


2v 2 COS &2 


4u 2 cos 62 


6w 2 cos 82 


{Jl cos 2 3 + J2 sin 2 3 ) 
* sin 2 83 + J3 cos 2 02 


(Ja - Ji) sin 2 
* sin 03 cos 03 


J3 COS 2 2 


e 2 











(Ja - Jl) Sill 02 

* sin 03 cos 03 


Jl sin 2 03 + J2 cos 2 3 





e 3 


2v 2 


4u 2 


6w 2 


J 3 COS 2 2 





J3 



TABLE VI: The matrix i^'J^W sin 6 2 (see text). 



Pv 



Px 



Li 


L 2 




















cos 0^ 


sin 9 3 


Ji 


Jl 


sin 02 sin $3 


sin 09 cos #3 




Ji 


cos f?2 cos #3 


cos #2 s i n ^3 




J 2 



9 
_0_ 

Ox 
_d_ 

00 



sill #2 1 s ' n 

— jtf 4 H 



2v 2 ' 4u 2 
2 sin 0o 



2 sin fl 2 



Sin f^2 

2u 2 



+ 



4 sin 02 
4^ 



3 sin 02 
6 sin 02 



+ 



_ sin 2 
2 sin 9 2 



a 

ee 3 




TABLE VII: The lowest 6x6 sub-matrix of the inverse of the 
matrix Q given in the Table EI The upper part is diagonal in 
the present approximation. 





p<p 


Px 


Pf> 


Li 


L 2 


£3 


Pv 


2v? ' 4uJ 


2 


3 
4^0 








1 

4"o 


Px 


2 


^T + A 

2u 2 4ug 


6 








2 


p« 


3 


6 
4~"o 


-±T + A 

2u> 2 4u 2 








3 

4^1 


Li 











1 

Jl 








u 














1 

Ja 





u 


1 


2 


3 
4^ 








1 

4^0 



with the matrix 1/ given by the eq. I|C4(1 . It is also con- 
venient to define the quantum operators 



Px 
Po 



—id I dip 
-id I dx 
-id I d-d 



(C7) 



to obtain 



d d d d 



d 

d8 2 



d_ 

de 3 



dip' d\' 01?' d&\ 
W (p v , p x , p#, L u L 2 , £3} 



(C8) 



The introduction of a set of conjugate momenta pk 
such that 



•JL 



kPi 



(C9) 



in the eq. I|C2|I gives 



+iG -^ dG1/ Y WVv ** 



(CIO) 



In our case, .4 = , &i = WQ{W, so that 

waw = w (w^g^w- 1 ) w = eq 1 (en) 



and the first term of the eq. (|C10fl takes the form given 
in the Tabic IVIII The second term of the eq. IICIOI) 
vanishes. In fact, only the derivatives with respect to 
82 or 83 could give a contribution to the sum, as the 
other variables do not appear in the elements of the sub- 
matrix. Moreover, the determinant G\ is simply propor- 
tional to sin 2 82 (with the proportionality factor depend- 
ing on dynamical variables which are outside the present 
sub-space), and the term G 1 / 2 in the equation l|C10|) 
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can be replaced by sin 6*2- The matrix i{G' 1 )~ 1 W sin 02 
is given in the Table IVll The sum of the derivatives of 
each element of the 5th row, with respect to 02 , plus the 
corresponding one of the 6th row, with respect to #3, 
would give the coefficient of the corresponding momen- 
tum operator, but it is easy to verify they cancel each 
other. 



APPENDIX D: THE POSSIBLE FORMS OF THE 
POTENTIAL TERM 

The form of the approximate potential-energy expres- 
sion depends on the details of the underlying microscopic 
structure that the model should try to simulate. There 
are, however, some general rules to which the expression 
of the potential energy must conform: it must be invari- 
ant under space rotation, time reversal, and parity. 

Irreducible tensors of different rank have been con- 
structed with each of the basic tensors a^ 2 ' and a^- 3 \ 
and scalar products of tensors of equal rank have been 
considered. In order to be invariant under time rever- 
sal, each term must contain an even number of octupole 
(3) 

amplitudes . Moreover, due to symmetry, only ten- 
sors T^ T p — [a( A ) ® a^]^ of even rank K can be ob- 
tained with the coupling of two identical tensors a^ x \ 
There are therefore two independent invariants of order 2, 
( a ( 2 )- a ( 2 )) and (a^-a^), two of order 3, (a^ -t£ } ) and 
( a ( 2 > • T33') , and seven of order 4, namely (T 2 ( f ] ■ t!£ ] ) 

(K = 2,4), (4P ■ 4P) (K = 2, 4, 6), and . 1^) 

(K = 2,4). Additional fourth-order invariants of the 
form ([«( 2 ) $ a ( 3 )]W • [a^ ® a^}^) (K = 1 - 5), are 
not independent from the above ones, and can be ex- 
pressed as linear combination of them with the standard 
rules of angular momentum recoupling. 

Invariant expressions up to the fourth order in the am- 
plitudes (A = 2,3) are shown in the Table IVTTT1 in 



terms of the dynamical variables /?2, 72, 73, w 2 = 
(A" 2 + Y 2 ), £, and v 2 = {ry 2 + ( 2 ). Expressions corre- 
sponding to different choices of the dynamical variables 
can be easily obtained. 

Here, only the terms up to the second order in the se- 
ries expansion of the "small" (non-axial) amplitudes are 
given. In this approximation, the fourth-order invariants 
built with the quadrupole amplitudes (C±, C2) result to 
be proportional to each other and to the square of the 
corresponding second-order invariant (A 2 ). Moreover, as 
it is clear from the Table IVTTT1 the three fourth-order in- 
variants built with the octupole amplitudes (C3, C4, C5) 
and the square of the corresponding second-order invari- 
ant (A?,) are not linearly independent of one another, and 
provide only two independent relations. 

One can finally observe that the variables X, Y, tj and 
£ always appear only in the combinations (A 2 + Y 2 ), 
(ry 2 + C 2 )- As long as the expression of the potential en- 
ergy does only depend on the invariants up to the fourth 
order, the angles i9 and p defined in the Eq. I|17|) are 
not subject to a restoring force. The situation is more 
complicated for the three variables related to the yi = ±2 
components of a^. Also in this case, however, it is possi- 
ble to construct combinations of invariants that contain 
only linear combinations of the squares of the variables 
f, 7 = \/572 - 73, and 70 = /?f 72 + V5f3 2 -f 3 - Moreover, 
in this case, the first two of them do appear in the com- 
bination /32/337 2 + 2(/3| + 5/3f)£ 2 . It is therefore possible 
to imaging a situation in which the potential energy is 
independent also of the angle x (although this is is not a 
direct consequence of the model, as it is the case for the 
angles $ and tp). 
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TABLE VIII: Invariants up to the fourth order built with the tensors and the limit close to axial symmetry. Here 

= [a K CS> a K ] ' L \ and we use the definitions v 2 = rj 2 + £ 2 , w 2 = X 2 + Y 2 . "Small" (non-axial) terms are approximated 
up to the second order. 

"gj „ c 2 ^ .. v 2 
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/3 2 

^2 
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B 2 
P3 






+ 10io 2 
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4 V 22 / 
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2 t J 2 12 
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"4/3 3 4 7 3 2 
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pIpI 
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+ 1373 2 ) 


+58iw 2 
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,}o04 2 i -i^L v 2 
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